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We report on recent advances on the computation of the matrix elements of the electroweak penguins Q7 and Qs 
which are relevant for the A J = 3/2 contribution to in the chiral limit. The matching of scale and scheme 
dependences between Wilson coefflcients and these matrix elements is done analytically at NLO in as- 
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where iAj = a^j exp(iSj) are K —> tttt isospin in- 
variant amplitudes and Si are final state inter- 
action (FSI) phases. The ratio Re(ao)/Re(a2) — 
22.2 is an experimentally well known quantity. 

In the limit m„ = md and ol^ed ~ ^'^'^ 
neglecting the tiny electroweak corrections to 
Re(a2), one gets 



Im(a2) 
Re(a2) 
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Im(e^GE) 
G27 



(2) 



including FSI to all orders in CHPT and up to 
Oip"^) non-FSI corrections ||l|,||J|]. The coupling 
G27 modulates the 27-plet operator describing 
K ^ -KIT at 0{p'^) in CHPT. Its value can be 
obtained from a fit of if ^ tttt and K —f tttttt to 
0(p4) amplitudes [§. 

The coupHng Ge appears in CHPT to 0{e'^p°) 
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CF^e^GE (C/^QC/)23+0(/)(G27, • • ^(S) 



The constant C was defined in [Q , _Fb is the pion 
decay constant in the chiral limit, Q is a 3 x 3 
matrix collecting the three light quarks electric 
charges and U = exp(zV2$/Fo) with $ a 3 x 



3 matrix collecting the octet pion, kaon, and eta 
pseudo-scalar boson fields. 

In the Standard Model, there are just two op- 
erators contributing to Im(e^GE); namely, the so- 
called electroweak penguins, Q7 and Qg. In the 
chiral limit, these operators form a closed system 
under QCD corrections. Its anomalous dimen- 
sions mixing matrix is known to NLO in the NDR 
and HV schemes [^,0. 

There has been recently a lot of work de- 
voted to calculate Im(e^GE), both analytically 
§,|,|l^,0|l| and using lattice QCD [pl |il,[l5| ,|l6| . 
Here, we report on the work presented in [ 1^ and 
to which we refer for explicit formulas and further 
references. For a review of the method we used 
and the matching procedure see [Q. Using this 
method, one can write exact results for the cou- 
pling Im(e^GE) in terms of integrals of full two- 
point functions in the chiral limit. They can be 
related to spectral functions via dispersion rela- 
tions and resummation of the effect of all higher 
dimensional operators in the OPE of the relevant 
two-point functions. 

The coupling Im(e^GE) can be written as 

- ^Fo^Me^Gfi) = -GlmCrim) {0\Q7\0)^{^^R) 

+lmCsim){0\Qs\0)xil^B)- (4) 

In the chiral limit and in the NDR scheme (for 
the HV scheme expressions see ll^), the matrix 
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elements above are 



with 
and 



" dU^ Inf^) -ImnJiiW (7) 
VMi?7 1" 



So and sq are the onsets of local QCD dual- 
ity. These results are exact in the chiral limit. 
The two-point function ^^j^{Q^) is the trans- 
verse part of (0 1 (F'' - ) (2;) 4- ) (3) (0) 1 0) 
in momentum space, and the two-point function 
^^ss+pp(Q^) is the singlet minus triplet combina- 
tion of {0\{S + iP){x){S-iP){0)\0) in momentum 
space, both in the chiral limit. Sec jl7j for their 
explicit expressions. In this reference, we show 
how scheme and scale dependences of the Wilson 
coefficients match analytically the ones of the ma- 
trix elements. This was done there explicitly in 
the NDR and HV schemes at 0(q;|). 

The first term in comes from a disconnected 
diagram and is order . The second term comes 
from connected diagrams and is order N^. We 

also know that ^'ss+ppiQ'^) OPE which 

starts at modulated by 0(a|) coefficients 

times (OlQvlO) and (OlQslO). As a consequence, 
it fulfills exactly a first Weinberg sum rule-like 
and very approximately a second one [l7| . 
Apart of fulfilling the two Weinberg sum rules- 
like, notice that the kernel of the sum rule (||) we 
are interested in has zeroes at t=0 and at t = 
fj.'^ which we chose to be = Sq- One expects 
therefore that keeping just the first pole in each 
channel is a good approximation for estimating 



the leading OPE behavior as happens in 7r+ — 7r° 
electromagnetic mass difference. 

In Q we discussed this sum rule using the 
known pion and the rji poles and including the 
first pion prime as a narrow width. The tt' con- 
tribution can be improved using a Breit-Wigner 
shape and the results do not change much due 
to the phenomenologically small coupling of the 
tt' . The scalar counterpart is more delicate and 
more work is needed. But as a first estimate, 
we used two hadronic models that ful- 

fill the QCD short-distance constraints and pro- 
duce values for L4, Lq and compatible with 
phenomenology. The scalar form factor in ||2^ 
is obtained from data and dispersion relations 
up to 1 GeV and Breit-Wigner shapes above. 
The result of using these models agreed with 
the results of naive narrow widths for the lowest 
scalar resonances. These were constructed to ful- 
fill the short-distance QCD constraints and also 
produced reasonable values for Lq and £§• Here, 
we have also tried Breit-Wigner shapes for the 
scalar mesons instead of narrow widths and again 
find results in the same ball park. Now, we also 
have used the scalar form factors obtained in ||2^ 
where the lowest scalar triplet and singlet reso- 
nances are generated dynamically for energies up 
to 1 GeV and Breit-Wigner shapes above and we 
get similar results. 

In all the estimates, we got negative correc- 
tions to the first term in (H) in the region between 
— 10% to —30%. Though the scalar-pseudoscalar 
sum rule (|^) cannot obviously be used at a quanti- 
tative level at present, the results above indicate 
that is difficult to have corrections larger than 
±30% to the first term in (||). 

The chiral hmit OPE of UIr{Q'^) starts at 
modulated by known coefficients times 
(OlQyiO) and (OlQslO). Using this, one arrives at 

§,raniiT7i 
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{0\Qs\0)^''''{so) . (9) 



This is the relation that we use for the quantita- 
tive analysis of (OlQslO). 

This sum rule and the one in can be calcu- 
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lated using the excellent ALEPH (23| and OPAL 
24 1 tau data. For the details of how we use the 



tau data we refer to |17 . Only to mention here 



that we have generated around 100,000 tau data 
distributions according to their covariance matrix 
and assigned to each one a single value of so- This 
value is the highest one allowed by data where 
the second Weinberg sum rule (WSR) is fulfilled. 
Both sum rules and ^ have as kernel. This 
makes the region above 1.5 GeV^ of increasing 
importance with respect to the low energy re- 
gion. We have chosen the second WSR since it 
differs from the one we are interested in just by 
one power less in the kernel. In addition, the sum 
rule in (|^) has a zero at t = fj.'^ = sq. We also used 
the highest value available by the data -which is 
always between 2 GeV^ and 3 GeV^- since one 
also expects QCD-Hadron duality to work better 
there. 

Combining the ALEPH and OPAL results, we 
obtain 



M2 = -[1.9 ± LO] • 10^3 GeV*^ 



(10) 



and always A/3 > in agreement with the recent 
analysis in [|2|. 



If we use Imll 
(01 gg| 0)^(2 GeV) 
|25[, we get 



(0-3) 
SS+PP 



it) 



in (§) and 



-(0.018±0.004) GeV^ from 



M2 = -[2.0 ± 0.9] • 10"^ GeV^ 



(11) 



This result is very compatible with the FESR 
analysis result we got in (p^). 

For our final result in the NDR scheme, using 
as(2 GeV) = 0.32 (see more details in [0), we 
quote 

(0|g8|0)^™(2GeV) = 
(1.20 ± 0.60 ± 0.15) • 10"^ GeV^ = 

(1.2±0.7) • lO^^GeV^ (12) 

Where the first error is purely experimentally and 
takes into account both ALEPH and OPAL re- 
sults as well as a possible variation of the local 
duality onset So and the second one is from the 
unknown 0{a^) terms in (|^) assuming a geomet- 
rical series. 

The sum rule in ^ is much better behaved 
and with smaller error bars due to to the zero at 



t ^ — sq. Using again the same strategy we 
explained above for the analysis of sum rule (^, 
we get 



y^LJ?(2GeV) = (4.35 ± 0.50) • lO"'' GeV 



(13) 



combining ALEPH and OPAL results. See (Tt) 
for further details. 

Gomparison with other recent determinations 
is made in Table |l|. For the results in 
the cases B^(2GeV) = B^(2GeV) = 1 and 



ImH 



(0-3) 



-ss+pp(t) = we used {0\qq\0)jjg (2 GeV) 
= -(0.018 ± 0.004) GeV^ from (2|, which is in 
agreement with the most recent sum rule deter- 
minations of this condensate and of light quark 
masses -see for instance- and the lattice light 
quark masses world average in [ |27[ | . 

Within the present accuracy of (Ol^glO), the dis- 
connected contribution to (^) -third line in Table 
|]- is perfectly compatible with our full result - 
fourth line in Table |l|- as well as the results from 
[|9|, p^ , p^ , p^|JT5[| , so that we cannot conclude a 
large deviation from the large Nc result within the 
present accuracy. Notice that we include in this 
result -third line in Table |l|- 0{as) corrections 
that are indeed leading order in 1/A^c (see (||) 
and [|l^) which are usually disregarded in the fac- 
torization approaches, this makes the chiral limit 
i?g (2GeV) parameter larger than one by around 
20% to 30%. 
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